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and b — ► dj, concentrating on non-standard models with new charged Higgs 
interactions. We find that the charged Higgs contribution to the asymmetry 
^ . for b — > S7 is small in such models due to a universal cancellation mechanism. 

The asymmetry is therefore difficult to distinguish experimentally from the 
Standard Model (SM) value, which is also small. The cancellation mechanism 
is found to be rendered inoperative in supersymmetry due to the presence of 
chargino loops. Unlike b — > sj, the rate asymmetry for b — > d*y in Higgs models 
can be quite different from its SM value, generally ranging from —20% to 
+20%. Specific model calculations are performed for the Three-Higgs Doublet 
Model and the "Top" Two-Higgs Doublet Model to serve as illustrations. We 
also offer some suggestions that may be helpful to experimentalists in the 
detection of the inclusive mode b — > dj. 
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I. INTRODUCTION 



The e + e~-based 5-factories are all performing quite well. Each of them should soon be 
producing of order 10 7 B-B pairs per year in a very clean environment. In another few 
years it is likely that one or more of these machines will yield of order 10 8 B-B pairs per 
year. 5-experiments at the hadron machines - Tevatron (BTEV), Hera-B and LHC-B - 
could increase this number by another 1-3 orders of magnitude. Therefore, studies of rare 
B-decays will continue to intensify for the next several years. 

Radiative B decays have been the subject of much theoretical and experimental interest 
over the past several years. While b — > 57 is loop-suppressed within the SM, it nevertheless 
has a relatively large branching ratio, allowing for a very important comparison between 
theory and experiment. A recent measurement of the branching ratio by CLEO PJ yielded 
the value B(b -> 57) = (3.15±0.35±0.32±0.26) x 1(T 4 , in good agreement with the ALEPH 
measurement and the SM prediction |^|-§J. The good agreement between the experimental 
and theoretical branching ratios places a strong constraint on many non-standard models. 
The recent CLEO measurements [p]] have also placed a rather modest bound on the partial 
rate asymmetry (PRA) in b — > 57, yielding -9% < A^p 1 < 42% at the 90% confidence 
level. This bound does not yet encroach on the very small PRA (of order 0.6%) predicted 
by the SM 0. Improved measurements of A h ^p sl are expected over the next several years 
at the 5-factories. These measurements should provide a powerful probe of new physics, 
particularly of models that contain non-standard CP-odd phases. 

A mode closely related to b — > 57 that is also of great interest is b — > dj. The branching 
fraction for this mode is expected to be of order 10~ 5 , so a good measurement of the branch- 
ing fraction should be possible with the 10 8 B-B pairs per year expected at the 5-factories, 
provided the difficulties in the experimental detection of b — > d'y can be overcome (see below). 
The PRA for b — > d'y is of order —16% within the SM and may also be observable within the 
next several years. It is important to note that the SM PRA for 6 — >■ g?t / is statistically more 
accessible than that for b — > 57. That is, the increase in the PRA for b — > d'y relative to that 
for b — > S7 more than compensates the suppression of its branching ratio. This statement 
is relatively robust and holds true even after making reasonable concessions for detection 
efficiencies in the two cases, as we discuss in Sec. [V (1 The statistical accessibility of A b ^p dl 
is quite intriguing and leads us to a serious consideration of this observable in the present 
paper. A discussion of some practical experimental techniques for distinguishing b — > d'y 
from b — > s'y may be found in Sec. 0. 

Our goal in the present work is to assess the sensitivity of the rate asymmetries noted 
above to effects coming from various non-standard models. Our main emphasis is on effects 
due to the exchange of charged Higgs bosons, but we offer insights regarding other models 
as well. A summary of our main results as well as those of some previous authors may be 
found in Table |. Previous studies of the b — > S7 PRA due to charged Higgs exchange have 
generically found quite small values fl7HT0"f . In this paper we will show that this trend is 



due to a partial cancellation between terms in a general expression for the rate asymmetry. 
The cancellation is a nearly universal feature of models containing only new charged Higgs 
interactions, and typically constrains b — > 57 rate asymmetries in such models to be at or 
below the few-percent level. Experimentally, such values would be difficult to distinguish 
from the (similarly small) value expected within the SM. A contrasting situation is found 
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for b — > d'y ||. In this case the PRA can differ appreciably from the SM result, in general 
ranging between —20% and +20%. 

The outline of the remainder of the paper is as follows. In Sec. [0] we review the general 
formalism used to calculate the PRA for inclusive radiative b decays. In Sec. |H] we focus 
specifically on Higgs models and show that A b ^p sl is destined to be quite small in such 
models, while A b ^p dl can be both large and different from its SM value. These general 
features are illustrated by calculations in a few specific Higgs models, namely in the Three- 
Higgs Doublet Model (3HDM) with Natural Flavour Conservation and the Top-Two-Higgs 
Doublet Model (T2HDM). Section |TV| contains a general discussion of supersymmetric effects 
on the asymmetries and shows how supersymmetric models manage to evade the accidental 
cancellation that Higgs models face. Section [V| contains a discussion of some practical 
experimental techniques that could be useful in distinguishing between b — > d'y and b — > sj. 
In Sec. |VT] we offer a few concluding remarks, as well as a summary of our results. 



II. GENERAL FORMALISM 

Let us begin by reviewing the general formalism used to calculate the PRA in inclusive 
radiative B decays. Further details of the calculation may be found in Ref. . The PRA 
in the inclusive decay b — > sj is defined by the ratio 

A ^ S1 = r(E^x sl )-r(B^Xr?) m 

CP T(B^X sl ) + T(B^Xrf) ' 

The next-to-leading order expression for the asymmetry may be found in Ref. |9j and is 
expressed in terms of the Wilson coefficients Ci as 

4p S7 = {^Im [C 2 C* 7 ] - ^ [v{z) + b(z, 6,)] Im [(1 + e s )C 2 C* 7 ] 

~Im [C S C* 7 ] + ^b(z, 5 7 )Im [(1 + e s )C 2 C* 8 ]j , (2) 
where z = m 2 c /m1 ~ (0.29) 2 and 

6, = ^--A 2 (p-^). (3) 

Vts V tb 

The second term in Eq. ([3]) gives an approximation to e s , written in the usual Wolfenstein 
parameterization, with A ~ 0.22. The functions v(z) and b(z, <5 7 ) are defined in Ref. 0; <5 7 
is related to the experimental cut on the photon energy: E 1 > (1 — <5 7 )i? 7 nax . An analogous 
CP asymmetry can be constructed for the decay b — > d'y and gives an expression similar to 
Eq. (H), but with the replacement e s e^, where [§] 

Kagan and Neubert have evaluated the coefficients of the various terms in Eq. (HP for dif- 
ferent values of the photon energy cut-off, 5 T They have also performed a more elaborate 
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calculation to take into account the effects of Fermi motion . For our purposes it is suf- 
ficient to neglect the Fermi- motion effects and to choose a particular value of S T Choosing 
<5 7 = 0.30, and taking a s {m b ) ~ 0.214, we find 



.4 



b— >s(d)7 
CP 



1(T 



ICVI 2 



1.17 x Im [C 2 C;] - 9.51 x Im [C 8 C 7 *] + 0.12 x Im [C 2 C*] 



-9.40 x Im 



e s{d )C 2 (C* 7 - 0.013 CI) 



(5) 



The rather large coefficient in the second term has led to the observation that models with 
enhanced chromomagnetic dipole operators could give rise to significant changes in the rate 
asymmetries, compared to those of the SM || . 

The Wilson coefficients in Eq. (||) are determined using the usual Renormalization Group 
Evolution (RGE) procedure PI and are to be evaluated at the h quark mass scale. The 



(6) 



results at leading order are JT2 

8 

C 2 {m h ) =J2 k *V ai 



C 7 (m b ) = v W/23 C 7 (mw) + f (r/ 14/23 - ?? 16/23 ) C 8 (m w ) + C 2 (m w ) £ hfl* (7) 

d i=l 



C 8 {m b ) = fj u / 23 C 8 {m w ) + C 2 {m w ) £ h^, 



(8) 



where r] = a s (m^)/a s (mj,). The coefficients fc 2 i, ^ and /ij are tabulated in Ref. |ll|]. The 
SM expressions for the Wilson coefficients are given by 



Ci M (m w ) = 1 



Cf\m w ) = —A(xt) 



Ci M (m w ) 



-D{x t ), 



(9) 
(10) 

(11) 



where x t = m^/rrv^ and where A and D are the well-known "Inami-Lim" functions [13]: 



A(x) = x 
D(x) = x 



3x 2 + 5x — 7 (3x 2 — 2x) \nx 



12(x - l) 3 
x 2 — 5x — 2 



+ 



2(x - l) 4 
3x In a; 



4(x-l) 3 2(x-l) 



(12) 
(13) 



Deviations from the SM may be incorporated into the Wilson coefficients at the W mass 
scale and can then be run down to the b mass scale using the RGE equations given in Eqs. @- 
(|). This procedure yields the following simple expressions for C 7 {m b ) and Cg(m b ) f|, 



C 7 {m b ) ~ -0.31 + 0.67 C 7 new (mvK) + 0.09 C™ w (m w ) 
C 8 (m b ) ~ -0.15 + 0.70 C 8 new (m w ), 



(14) 
(15) 
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where we have used fj ~ 0.56 (corresponding to a s imz) = 0.118) and where the Wilson 
coefficients at the W mass scale are given by 



C 7 (m w ) = C* M (m w ) + C*™(m w ) 
Cz(m w ) = Cf A (m w ) + Cl 



lm w 



(16) 
(17) 



Inserting Eqs. ( |H| ) and ( |T5| ) into Eq. (||) yields the following general expression for the CP 
asymmetry 



.4 



b— >s(c()7 
CP 



10" 



\C 7 (m b )\i 

+3.21 x Im 



{-1.82 x Im [C? w ] + 1.72 x Im [C 8 ncw ] - 4.46 x Im [Q cw Cf w *] 



s(d) 



2.18 



0.26 C 8 new *) 



(18) 



Note that if C^ 6 ™ and Cg ew are approximately equal in a particular scenario, then the third 
term in the above expression is close to zero, and the first two terms nearly cancel one 
another. The main contribution to the PRA in that case is the last term. Such a scenario 
yields a rather small asymmetry for b — > 57 (since the asymmetry is essentially proportional 



CS 



111 



to e s ), but possibly a sizable asymmetry for b — > dj. As we shall see below, Cy ei 
multi-Higgs doublet models, so that b — > 57 generically has a small rate asymmetry in such 
models. 

Setting both Cj cw (m\y) and Cg ew (mw) to zero in the above expression yields the SM 
predictions for the CP-violating rate asymmetries in b — > S7 and b — > d'y, 



A 



b— >s(d)7 
CP 



0.334 x Im 



or, in terms of the angles of the familiar unitarity triangle |14 

-0.334 x 



/,-.-. _ 334 x y2 sin Pckm sin 7 Ab ^ dj _ n nn A sin a sinfe M 



A CP 



sin a 



4 ~ 
A CP 



sm7 



(19) 



(20) 



Substituting the approximate expressions for e s and from Eqs. @ and @ and using the 
"best-fit" values (p, 77) = (0.20,0.37) |T§, we find 



A CP 



SM 



0.( 



A 



b— >d7 
CP 



SM 



-16% 



in agreement with the estimates of previous authors 0]. The rate asymmetry for b — > 57 
is potentially a sensitive probe of non-standard physics, since its value within the SM is so 
small. Finding an experimental value larger than a few percent would be a strong indication 
of new physics. For b — > d'y, the rate asymmetry is already large within the SM; in this case, 
contributions due to non-standard physics might be more difficult to disentangle from the 
SM contribution unless they lead to an appreciable change in the PRA. 



III. MULTI-HIGGS DOUBLET MODELS 

One of the simplest ways to extend the SM is to expand its Higgs sector by adding 
extra Higgs doublets. Despite their apparent simplicity, Multi-Higgs doublet models can 
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have rich phenomenological consequences. Of particular interest within the context of the 
present work are those models that give rise to new sources of CP violation. If one imposes 
Natural Flavour Conservation (NFC) ||16|| , CP violation in the charged Higgs sector does 
not appear until two extra Higgs doublets have been added. The Three-Higgs Doublet 
Model (3HDM) with NFC, incorporating spontaneous and/or explicit CP violation, has 
been studied extensively over the last several decades [pT7|— |2lf] and in many instances it 
leads to large CP-violating effects. Another option is to relax the requirement of NFC. In 
this case, CP violation appears with only one extra Higgs doublet. The so-called "Type- 
Ill" Two-Higgs doublet models can lead to interesting phenomenological consequences . 
One variant on this class of models - the "Top-Two-Higgs Doublet Model" (T2HDM) - 
gives special status to the top quark p3| , p4| . The T2HDM has couplings unlike most other 
Higgs models and can produce significant effects on CP-odd observables in decays such as 
B — > ipK s and B ± — > ipK ± and in D-D mixing |f23l-p5f . 

Bearing in mind the above comments, it would seem natural to assume that Multi-Higgs 
doublet models could give rise to significant effects in the rate asymmetries for radiative h 
decays. This intuitive guess is correct for b — > dj, but is not correct for b —>■ 37. Small 
values for A^p 1 were reported in previous studies f7]-|T0l|, but the reason for the smallness 
of the asymmetry was not elucidated. In the following we show that the rate asymmetry for 
b — > s'j is typically less than a few percent in Multi-Higgs doublet models. The analogous 
asymmetry for b — > d'j can nevertheless be strongly affected by the presence of extra Higgs 
doublets. Our general arguments will be supported by model calculations in the 3HDM and 
T2HDM. 

The rate asymmetry for b — > 57 is small in Higgs models because 
Cg ew (my/) I 'C" ew (m w )r^l in these models. As noted in the discussion following Eq. (|18|). 
this approximate equality leads to a strong suppression of the rate asymmetry. If we only 
consider contributions to the rate asymmetry arising from diagrams containing a charged 
Higgs and a top quark in the loop 



, the "new physics" corrections to C7 and C% have the 



following general form 



Cr»(m w ) 



f A A(y t ) + f B B(y t ) 
f A D(y t ) + f B E(y t ), 



(21) 
(22) 



where y t = mf/m 2 H , and A(y) and D(y) are the Inami-Lim functions introduced above. The 
functions B(y) and E(y) are given by 



B(y) 
E(y) 



y 



% - 3 _ (3y - 2) lay 
\ny 



12(2/ - l) 2 

y-3 



+ 



4(y-l)2 2(j/ -1) 



(23) 
(24) 



Figure |l|(a) shows the various Inami-Lim functions, plotted as functions of the dimensionless 
variable yt, while Fig. |l|(b) shows plots of the ratios Dj A and E/B. As is evident from these 
plots, D/A and E/B are both of order unity for Higgs masses above 100 GeV. As a result, 
Cg ew /Cy ew is of order unity and the CP asymmetry for b — > 37 is expected to be rather 
small. Note that the deviation from unity is slightly larger for D/A than for B/E when 
mil ~ Tat- Thus, slightly larger asymmetries may be expected in scenarios where Ja^Jb- 
(See the footnote and discussion following Eq. (|||).) 
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Let us quantify the above argument by parametrizing the deviation of Cg ew /C^ ew from 
unity in terms of a complex quantity (. Defining 



Cr(mw) = (1 + C)C 7 ncw Ky) , (25) 



we have 



A b cp {d)1 - , ' (-0.10 Im [C 7 ncw ] + 1.72 Im [CC 7 ncw ] - 4.46 |C 7 ncw | 2 Im [(} 
6 7 m 6 2 I 



+3.21 Im 



1- 2.44 C™ w * - 0.26 CCT**) >■ (26) 



Each of the first three terms in the above expression is small, either because it is suppressed 
by ( or, in the case of the first term, because it has suffered an accidental cancellation. 
Since e s is itself quite small, the asymmetry for b — > sj will never be much larger than 
several percent. The present analysis is quite general and predicts that the CP asymmetry 
in b — > S7 will be quite small in Higgs models. This result is independent of whether or not 
the corrections to the Wilson coefficients have complex phases. The situation for b —>■ d^f is 
quite different. While it is still true that the first three terms in Eq. ( |26| ) make only a small 
contribution to the asymmetry, the last term now contains e^, which is relatively large and 
can have a significant effect on the asymmetry. 



A. Three-Higgs doublet model with NFC 



We now consider the 3HDM with explicit CP violation, since the 3HDM with spontaneous 
CP violation alone is disfavoured by data |T^^T| • The charged-Higgs Yukawa couplings for 
this model may be written as flTj 



(2V2G F )^H + (XU L VM D D R + YU R M V VD L + ZN L M E E R ) + H.c. 



(27) 



where X, Y, and Z are complex parameters |[27|| . Note that the substitution X, Z — > tan/3 
and Y — > cot f3 yields the Yukawa interactions of the Type-II 2HDM . The Wilson coefficients 
at the my/ scale are given by 



C 7 (m w ) = --A(x t )-XY*B(y t ) 
Cz(m w ) = - l -D{x t )-XY*E(y t ) 



y\ 2 A{y t ) 
Y\ 2 D(y t 



(28) 



Thus 



f A 
SB 

c 



-V1 2 

6 1 1 
-XY* 



XY*[E(y t ) - B(y t )} + \\Y \ 2 [D(y t ) - A(y t )] 



(29) 



XY*B{y t ) + \\Y\>A{y t ) 
In the numerical work that follows, |£| is typically found to be less than 0.20, so Eq. (|26| 



yields a small asymmetry for b — * 57, as noted above [28 
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The parameter space for this 3HDM depends on four independent quantities, which may 
be taken to be \X\, \Y\, arg(Xy*) and m#. The phase of Y may be absorbed into the 
definition of X without affecting any observable in b — > s(d) r y. We thus take Y to be real 
and non-negative and write X as 



X = X R + iX Il 



(30) 



where Xr and Xj are the real and imaginary parts of X. Possible experimental constraints 
on X and Y come from neutral kaon decays, the neutron electric dipole moment, B-B mixing 
and the (CP-averaged) branching ratio for b —>■ 57 p0y^|l8[ , j29| . Another constraint may be 
obtained by requiring the Higgs sector of the theory to be perturbative ||30fl . Most of these 
considerations lead to fairly mild constraints on the parameter space of the 3HDM. Pertur- 
bativity requires \X\ < 130 |30| , while B-B mixing requires \Y\ < 1 for m# ~ 200 GeV. 
The neutron electric dipole moment places a very weak constraint on lm(XY*) pT| , |32| . A 
much stronger bound may be obtained from the branching ratio for b — > 57, which yields, 
for example, lm(XY*) < 3.5 for m# = 200 GeV (see Eq. (^) below, or Fig. [|). ex may also 
be used to constrain the parameter space of this model. If p and 77 take their SM central 
values, the constraint on Y is similar to that coming from B-B mixing. Long distance effects 



due to charged Higgs exchange may also be important for [pT8|j33[1 . Given the theoret- 
ical uncertainties associated with long-distance calculations, we ignore such effects in our 
analysis. 

The strongest bounds on X and Y come from the (CP-averaged) branching ratio for 
b — *■ 57. In order to minimize uncertainties due to the mass of the bottom quark, it is 
common to define the ratio 



R 



r(6 -»■ s 7 ) 



r(6 — > cev) 
Calculation of R at leading order gives |34] 

R ~ 0.0231 x |C 7 (m b )| 2 , 



(31) 



(32) 



where C-j{rn^) is given in Eq. (0). The region allowed by the b — >• 57 branching ratio is 
bounded by two concentric circles in the complex A-plane for given values of Y and m# . 
This may be seen by writing 



CV( m &) = — c — Xd, 
where c and d depend only on Y and m#: 



(33) 



c = 0.31 + - [0.67A{y t ) + 0.09D(y t )] \Y\ 2 
o 

d = [0.67B(y t ) + 0.09E(y t )] Y. 
The constraint due to the measured rate for b — > S7 then takes the general form 



R mm < 0.0231 x d 2 



x R + 1 ) +Xf 



< R T ' 



(34) 
(35) 



(36) 



S 



that is, the allowed region is bounded by two concentric circles in the complex X-plane. 
i? max and R min depend on the experimental value for the ratio R as well as on various 
experimental and theoretical uncertainties. The above expression may be used to find an 
upper bound on |Im(XF*)|. In terms of the geometrical picture presented above, the upper 
bound occurs at the "top" and "bottom" of the outer circle, yielding 



1 / Dmax 

|Im(Xr*)| < — rW— • (37) 

1 V ;I ~ 0.675 (y t ) + 0.09E{y t ) V 0.0231 V ; 

A search of the parameter space for this 3HDM yields results consistent with the general 



arguments given for Higgs models at the beginning of Sec. pj. The rate asymmetry for 
b — > 57 is found to be smaller than three or four percent, while that for b —>■ d'j can 
vary between ±20%. Figures 2 and 3 show plots that are representative of the results 
obtained. In each figure, contours of constant asymmetry for b — ■> s(d)7 are shown in the 
(Xr,Xi) plane for fixed values of Y and m#. Labels for the various contours are given 
as percentages. The size of the ring-shaped allowed region is determined by the maximum 
and minimum values allowed for the ratio R. The current experimental value for this ratio 
may be determined by using the weighted average of the CLEO Q and ALEPH values 
for the b — > 57 branching ratio, B(B — > X s 7)| ave ~ (3.14 ± 0.48) x 10~ 4 , as well as the 
branching ratio for b — > ceu, B(B — > X c ev) ~ 0.104 ± 0.009 fl3"5 |. This procedure yields the 
value R = (3.02 ± 0.53) x 10~ 3 . In addition to the experimental uncertainty for this ratio, 
there is another theoretical uncertainty. Since we are using LO expressions for the Wilson 
coefficients, we estimate the magnitude of this uncertainty to be 5R theorj = 0.7 x 10~ 3 . 
Combining the experimental and theoretical uncertainties in quadrature and doubling the 
result to get a "2a" uncertainty, we obtain R min = 0.0013 and i? max = 0.0048. These are 
the values used in the plots. 

Figure ^| compares the asymmetries obtained for two different values of the Higgs mass 
with Y held constant, while Fig. || keeps m# constant and shows the effect of changing Y. 
While the size of the allowed region exhibits a strong dependence on the values chosen for 
Y and m# (note the changes in the horizontal and vertical scales in the plots), the contours 
themselves look very similar from plot to plot. This is particularly true in the case of b — > d'j, 
whose asymmetry contours appear to be exact replicas of one another once an appropriate 
scaling and shifting has been performed in the (Xr, Xj) plane. In order to understand the 
scaling behaviour of A b c ^f' y , it is easiest to return to one of the original rate asymmetry 
expressions, Eq. (|j). The first three terms in Eq. (|5[) are identical for b — > 57 and b — > dj in 
this model. According to the discussion following Eq. (|26D , these terms combine to give a 
contribution to the asymmetry of at most a few percent. A good approximation for b — > d'j 
is thus obtained by neglecting these terms altogether (as well as the small contribution 
proportional to e^C^Cg) and writing 

b ^ dl -9.40 x IP' 2 

CP |C 7 (m h )| 2 l e dC 2 C 7 ] . (38) 

Changes in Y or m# lead to changes in the functions c and d defined in Eqs. (0) and (p5|) . 
It is always possible to perform a linear transformation on Xr and Xj that compensates 
for this change and leaves Cj(mi,) invariant. Since both the above approximate expression 
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for A 1 ^ 1 and the expression for the rate, Eq. fl32|) , depend only on C-j(rrib 



they too are 

invariant under this transformation [3~o ] . The variations in the contours for b —>■ sj may be 
attributed to the fact that C^rrib) - which has a larger effect on b —>■ 57 - is not invariant 
under the transformation. 

Figure |] shows a few representative correlation plots demonstrating the relationship 
between the asymmetries for b — > S7 and b — > d^y. It also shows the relationship between 
the asymmetry and branching ratio for b — > ^7. The branching ratios in these plots have 
been normalized to the theoretical expression for B(b — > 57) within the SM. The ratio of 
branching ratios has the simple form 



B. 



3HDM 



(6 - d 7 ) 



C 7 {m b )\ 



3HDM 



C7(^6)|SM 



(39) 



The points for these plots are taken to be uniformly distributed on square grids inside the 
respective allowed regions in the complex X-plane. This method of choosing points leads to 
the observed "textures" in these plots. 



B. Type-Ill Two-Higgs Doublet Models 

In Two-Higgs Doublet Models with natural flavour conservation |H| (as in types I and II 
P?|), there exist no new CP- violating phases in the charged-Higgs Yukawa couplings beyond 
that of the Cabibbo-Kobayashi-Maskawa (CKM) |38j matrix. Consequently, the charged 
Higgs contribution to the b — > s 7 rate asymmetry is both CKM- and QCD-suppressed 
(the latter by a factor a s /n) and the asymmetry is small, as in the SM. The b — > dj rate 
asymmetry is also similar to its SM value of approximately —16%. In Type III 2HDM's 
|22|| , where no adhoc discrete symmetry is imposed, one generally expects new CP- violating 



phases in the Yukawa couplings (as well as flavour- changing neutral Higgs interactions). As 
a result, the CKM and GIM suppressions present in the SM for b — > 57 are in general no 
longer operative, and a potentially large partial rate asymmetry in b — > 37 seems possible. 
This expectation is invalid due to the accidental cancellation mechanism discussed above. 

As a concrete example of a Type-Ill 2HDM, we consider the 2HDM for the top quark 
(T2HDM) P3]J2^ ]. This model is designed to accommodate the heaviness of the top quark 
by coupling it to a scalar-doublet with a large VEV. The other five quarks are coupled to 
the other scalar-doublet, whose VEV is much smaller. As a result, tan j3 is naturally large in 
this model. A distinctive feature of the T2HDM is its tan /3-enhanced charm quark Yukawa 
coupling, with the consequence that the charged Higgs sector can significantly affect the 
neutral kaon system and the b — ► 57 rate. 

Although there are two new CP phases in the T2HDM from the unitary rotation of the 
right-handed up-type quarks, only one of them contributes to radiative b decays. This is 
the phase 8 of the complex parameter £ = |£|e _i<s , measuring the cr-Ir mixing p4| |. As was 



done before, we take |£| = 1 throughout this paper. In contrast with the 3HDM, the virtual 
charm quark contribution in the T2HDM can occasionally be of the same magnitude as that 
of the top quark. The Wilson coefficient for b — > (q = s, d) is given by p4| 



i=c,t 



tan 2 $ + — — f S T 1/*) ftan 2 /5 + l 
miV* q V h q V 
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x (flfo) + l -A{ yi ) [-1 + -J— (SV) (cot 2 /3 + 1)1 } , (40) 

where k 19 = —VibV* q /(V t bV t *), yi = (mi/m H ) 2 , and V denotes the CKM matrix. The defini- 
tion of the matrix E may be found in p4| . The expression for Cg ew (mw) is obtained from 
Cj cw (m w ) by the substitutions A(yi) — > and f?(yi) — > £J(yi). 

We have studied the rate asymmetries over much of the parameter space in this model and 
have found the PRA for b — > 57 to be very small, consistent with our general observations 
above. For b — > dj, the PRA shows a strong dependence on the CKM phase (see Eq. fl35|)) 
but is barely affected by the new phase 5. Interestingly, the CKM phase can take a wide range 
of values in this model because the observed CP violation in the neutral kaon system may 
come partly or solely from the charged Higgs sector [[24j. As a result, the PRA for b — > d*y 
can be very different from the SM expectation, and can even be quite small. These features 
are illustrated in Figs. [| and |[ The shaded regions in these plots indicate the "allowed" 
regions, while the white regions are ruled out by constraints coming from ex, ArriK and the 
branching ratio for b — > 57. The constraints coming from ex and Amx were discussed in 
detail in Ref. \24\ for several scenarios and have been adopted here without change for similar 
scenarios. Figure [| illustrates the case of a SM-like scenario: the Wolfenstein parameters 
p and 7] take their SM best fit values, and the CP asymmetry for b — > d^ varies in the 
neighborhood of its SM prediction. Figure |5] shows a real-CKM scenario, where 77 = and 
p = 0.42. In this case the CP asymmetry for b — * dj is close to zero. Note that the recent 
CDF measurement of the CP asymmetry in B — > ipK$ disfavours a negative value for 



r], both within the SM and within the T2HDM p4| . We thus estimate the maximum CP 
asymmetry for b — > dj within the T2HDM to be of order +4%, as indicated in Fig. ^. 

It is worthwhile to compare the rate asymmetries in the 3HDM and the T2HDM. It is 
evident from Figs. that the rate asymmetry for b — > S7 is small in both models, as argued 
above on general grounds. In contrast, the rate asymmetry for b — * d'y in these models 
can differ appreciably from its SM value, —16%, taking on values between ±20% in the 
3HDM and between approximately —16% and +4% in the T2HDM. An important difference 
between the two models is that A 6 ^ 1 in the T2HDM is almost completely determined by 
the CKM phase, while the analogous asymmetry for the 3HDM exhibits a strong dependence 
on the charged Higgs sector CP-violating phase. 



IV. SUPERSYMMETRY 



Supersymmetric predictions for the b — ► 57 rate asymmetry are highly model-dependent 
|5|j9|j40|-^6l . Restrictive models, such as minimal supergravity, yield asymmetries of order 2% 
or less when the electron and neutron electric dipole moment constraints are taken into ac- 
count |40 |. Relaxing some of the assumptions of minimal supergravity leads to asymmetries 
as large as 10% f4"l|, while allowing for non-CKM-like intergenerational squark mixing can 



lead to asymmetries of order 15% [0,0 or larger if the gluino mass is significantly lighter 
than the squark masses 0. In this latter case, gluino-mediated diagrams become important 
in addition to the usual chargino- and charged Higgs-mediated diagrams. Allowing for the 
violation of i?-parity can lead to asymmetries of order 17% [IB]. A recent calculation of the 
b — » c?7 rate asymmetry in a fairly restrictive supersymmetric scenario yielded values in the 
ranges -(5 - 45)% and +(2 - 21)% [|7|. 
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The b — > S7 rate asymmetry within supersymmetric models can be appreciably larger 
than that in Higgs models or in the SM. It is instructive to see how supersymmetric models 
evade the accidental cancellation encountered by Higgs models. In order to understand this, 
let us consider the chargino-squark loop contributions that must typically be added to the W 
and charged Higgs contributions considered thus far. These diagrams give rise to additional 



terms in the Wilson coefficients [48], which we denote by 



C^mw) = foKl{y) + f D K 7 2 (y) , (41) 
C£(m w ) = f c K*(y) + f D Kl(y) . (42) 

Here fc and fo parameterize the couplings and y is the square of the ratio of the chargino 
and squark masses. The chargino loop for CV involves two new Inami-Lim functions, K\ = 
6A(y) + 5D(y) and K\ = 6B(y) + 5E(y), whereas for C$ the two functions are given by 
Kf = 9A(y) - GD(y) and K\ = 9B(y) - 6E(y). Since D/A ~ 1 and E/B ~ 1, we have 
K\jK\ ~ 3/11 and K\jK\ ~ 3/11. As a result, C^{mw) is generically a few times larger 
than C%{mw)i an d there is no longer a cancellation between the first two terms of Eq. 
Thus, a large asymmetry in supersymmetric theories is possible. 



V. COMMENTS ON EXPERIMENTAL TECHNIQUES AND SENSITIVITIES 

The experimental detection of the inclusive process B — > (i.e., b — *> dj) is very 
challenging. The main difficulty is that the smaller signal for b —>■ dj must be isolated from 
the much larger "background" of b — > sj. This requires ensuring that the signal events for 
b — » d'-f do not carry net strangeness. In the following we discuss two ideas that may be 
used to determine the branching ratio and partial rate asymmetry for B — > X^j- We also 
compare the experimental sensitivity to the PRA's for b —>■ s'j and b —>■ dj. 



A. Multiplicity cut 

One strategy that may be helpful in determining B(B — ► X^) is to use a multiplicity 
cut so that one deals with a semi-inclusive sample consisting of a maximum of n (say 5) 
mesons. Using such a measurement along with the corresponding one for B — > A s 7, one can 
deduce the total inclusive rate for B — > Xd'y to a very good approximation via: 

B(B -> X dl ) 
B(B -> X sl ) 

T(B -> 7 + 2vr) + T(B -> 7 + 3vr) + ■ ■ • + T(B -> 7 + mr) 
~ T(B -»• 7 + K + vr) + T(B -»• 7 + K + 2vr) + • • • + f(S -»• 7 + AT + (n - l)vr) 



(43) 



Note that each value of n for which a measurement is available leads to such an equation. 
Thus, in principle, the measurements could be repeated for different values of n to help 
reduce systematic uncertainties [49]. There are two key assumptions used in Eq. (|43|) . The 



first is that SU (3) holds; this assumption could introduce errors at the 10% level. The second 
assumption is that non-spectator contributions may be ignored - i.e., Eq. (|43| ) assumes that 
the spectator approximation holds precisely. This latter assumption should work very well 
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for neutral 5's. For charged 5's, the annihilation graph can make a contribution of order 
20% [p0 | , so that Eq. ( fl"3|) will receive somewhat larger corrections for nevertheless 
Eq. ( f4"3| ) provides an excellent basis for the determination of an approximate expression for 
B{B^X dl ). 



B. End point study via an energy cut 

Another possibility is to place a mass cut on the recoiling hadron so that mx < {itl-k + 

ttik) ~ 635 MeV. Correspondingly there is a minimum energy cut on the photon, e.g., 

E 1 > 2.60 GeV This would ensure that the decay products result from b — > and 
not from b — > sj. PRA determination from this would be quite straightforward: 

M = B\B ^ X dl ) - B\B ^ X^) 

where B' stands for the branching ratio with the energy cut on the photon. 



C. PRA in b -> dry 

Since the number of 5's needed for searching for a given asymmetry Acp scales approx- 
imately as (BAqp) -1 , the PRA in b — > d'y is statistically much more accessible than that in 
b — > 57. The point is that, based on the SM, the PRA in b — > dj is expected to be 20-40 
times larger than that in b — > sj. Not only does this offset the difference in their branching 
ratios, it may even overcome the decrease in detection efficiency for b — > d^y. 

As an illustration, suppose B(B -> X d -f) ~ B{B -> A,7)/20 ~ 1.5 x 10~ 5 . Then, to 
establish a 16% PRA to 3cr significance with a 20% detection efficiency and a 25% tagging 
efficiency, the needed number of B-B pairs is: 

Nl a = 4.7 x 10 8 (45) 

B 1.5 x 10^ 5 x 0.2 x 0.25 x 0.16 x 0.16 v ; 

In contrast, let us take A^p 1 ~ 0.6%, a detection efficiency of 70% and a 50% tagging 
efficiency. In this case the number of B-B pairs needed is: 

Nl a = 2.2 x 10 9 (46) 

B 3.2 x 10- 4 x 0.7 x 0.50 x 0.006 x 0.006 v ; 

Therefore, establishing a PRA in b — > d*y may be less difficult than in b — > 57. 



VI. CONCLUDING REMARKS 

The study of CP violation in B decays is one of the central themes at the various B- 
facilities. In particular, the standard model CKM paradigm of CP violation will be tested 
and new sources of CP violation will be searched for. In this work, we have shown that the 
charged Higgs loop contribution to the rate asymmetry in b — > 37 is small in all models due to 
an accidental cancellation in the loop integrals. More specifically, in charged Higgs models, 
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Cg ew ~ C" ew . As a result, the first two terms in Eq. ([Tj|) nearly cancel and the third term is 
suppressed, leading to a small asymmetry. Therefore, any experimental discovery of direct 
CP violation in b — > 57 at the 10% level will be due to new sources of CP violation beyond 
charged Higgs, for example supersymmetry. We have also explained why chargino loops in 
supersymmetry do not undergo this cancellation. This is because Cj(mw) is generically 
a few times larger than C$(mw)- Consequently, sizeable asymmetries are possible within 
sup er sy mmet ry. 

The direct CP asymmetry in b — » c/7 behaves quite differently from its counterpart in 
b — > S7. In Multi-Higgs doublet models, the asymmetry can be large and of opposite sign 
relative to the SM. In the T2HDM, in particular, CP violation is mainly determined by 
the CKM phase, and the asymmetry can show a strong deviation from the SM prediction 
if the observed CP violation in the neutral kaon system has a sizable component from the 
charged Higgs sector. These features make the the experimental measurement of the rate 
asymmetry in b — > d'y very worthwhile. Our main results are summarized in Table | along 
with representative results from various other models. 

We have also offered some suggestions that we hope would lead to the measurement of 
the inclusive decay B(B — > X^j), and eventually even the PRA in that mode. 
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TABLES 



Model 


A b ^ (%) 


A h ^ (%) 


SM 


0.6 


-16 


2HDM (Model II) 


~ 0.6 


16 


3HDM 


-3 to +3 


-20 to +20 


T2HDM 


~ to +0.6 


16 to +4 


Supergravity [g|]4l],|47| 


10 to +10 


-(5-45) and +(2-21) 


SUSY with squark mixing p2|j43|,p| 


15 to +15 




SUSY with i?-parity violation |46|] 


17 to +17 





TABLE I. CP asymmetries in various models. Quoted ranges are approximate; see the text for 
details. The blank entries represent quantities that have not, to our knowledge, been considered in 
the literature. The rate asymmetry for b — > 57 from SUSY with squark mixing could be larger than 
the quoted ±15% if the gluino mass is significantly lighter than the squark masses [9]. Note that the 
range quoted for the i?-parity violating case assumes m VT ~ 10 keV. The asymmetry is negligible 
if m„ T <C 10 keV, as is indicated by the SuperKamiokande atmospheric neutrino oscillation data 
[52]. 
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FIGURES 
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FIG. 1. (a) The Inami-Lim functions as defined in the text, (b) Ratios of the Inami-Lim 
functions. 
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FIG. 2. Contour plots of partial rate asymmetries in the decays b — > S7 and b — ► d"f in the 
3HDM with p = 0.20 and r/ = 0.37. The shaded ring-like regions are experimentally allowed. The 
contour lines show the CP asymmetries as percentages. 
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FIG. 4. Correlation between the rate asymmetries of b — > 57 and b — ► c?7 in the 3HDM. Also 
shown are the ranges of the ratio of branching ratios for the radiative decays. 
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FIG. 5. Contour plots of PRA's in b — > s~f and b — ► in the T2HDM - SM-like scenario: p 
and T] take their SM best fit values (p = 0.20 and rj = 0.37). The experimentally allowed regions 
are shaded, and the contours show the asymmetries as percentages. 
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